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^ ; Abstract 

^ \ The full ADHM-Nahm formalism is employed to find exact higher charge caloron so- 

^ ■ lutions with non-trivial holonomy, extended beyond the axially symmetric solutions 

found earlier. Particularly interesting is the case where the constituent monopoles, 

that make up these solutions, are not necessarily well-separated. This is worked out 

^ \ in detail for charge 2. We resolve the structure of the extended core, which was previ- 

O \ ously localized only through the singularity structure of the zero-mode density in the 

^ ■ far field limit. We also show that this singularity structure agrees exactly with the 

Tj- ■ abelian charge distribution as seen through the abelian component of the gauge field. 
O 

O 



As a by-product zero-mode densities for charge 2 magnetic monopoles are found. 



1 Introduction 

Oh! 
(D . 

' Calorons are instantons at finite temperature. For a long time the influence of a 
background Polyakov loop on the properties of these topological excitations has been 
^ . neglected. Solutions were constructed long ago PP and were studied in detail in the 
semi-classical approximation j2]. In all these studies the Polyakov loop at spatial in- 
finity (also called the holonomy) was trivial, i.e. an element of the center of the gauge 
group. That the influence of the background Polyakov loop on the topological excita- 
tions can be dramatic is particularly clear in the confined phase, where on average its 
trace vanishes. Caloron solutions in such backgrounds were constructed only relatively 
recently [SI E] and can be seen as composed of massive monopole constituents with 
their magnetic charges adding to zero. 

It was observed that the one-loop correction to the action for configurations with 
a non-trivial asymptotic value of the Polyakov loop gives rise to an infinite action 
barrier, which were therefore considered irrelevant However, the infinity simply 
arises due to the integration over the finite energy density induced by the perturbative 
fiuctuations in the background of a non-trivial Polyakov loop The proper setting 
would therefore rather be to calculate the non-perturbative contribution of calorons 
(with a given asymptotic value of the Polyakov loop) to this energy density, as was 
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first successfully implemented in supersymmetric theories jB], where the perturbative 
contribution vanishes. It has a minimum where the trace of the Polyakov loop vanishes, 
i.e. at maximal non-trivial holonomy. 

In a recent study at high temperatures, where one presumably can trust the semi- 
classical approximation, the non-perturbative contribution of these monopole con- 
stituents (also called dyons) was computed jZj. When added to the perturbative con- 
tribution with its minima at center elements, a local minimum develops where the 
trace of the Polyakov loop vanishes, deepening further for decreasing temperature. This 
gives support for a phase in which the center symmetry, broken in the high temperature 
phase, is restored and provides an indication that the monopole constituents are the 
relevant degrees of freedom for the confined phase. 

Also lattice studies, both using cooling js8, and chiral fermion zero-modes |9| as fil- 
ters, have now confirmed that monopole constituents do dynamically occur in the con- 
fined phase. A charge 1 caloron is seen for SU(?2) to consist of n constituent monopoles. 
In the deconfined phase, due to the fact that the average Polyakov loop becomes a 
center element, the caloron returns to the form known as the Harrington-Shepard so- 
lution PJ. The latter can also be interpreted as consisting of constituent monopoles, 
however, with n — 1 of them being massless. 

To be precise, for self-dual configurations in the background of non-trivial holonomy 
the masses of constituent monopoles are given by Sn'^uj/P, with Uj = /ij+i — /xj. The 
fii are related to the eigenvalues of the Polyakov loop at spatial infinity, 

Poo = lim Pexp( / Ao{t,x)dt) = g^exp{27ridiag{ni,iJ,2,...,Hn))g, (1) 

►oo JO 

(this expression assumes the periodic gauge A^(t, x) = A^(t + P,x)) where g is the 
gauge rotation used to diagonalize Poo and P is the period in the imaginary time 
direction, related to the inverse temperature. The eigenvalues exp(27ri/ij) are to be 
ordered on the circle such that /ii < yU2 < • • • < /^n < yUn+i, with /in+j = 1 + yUj and 
Y^7=if^i = O5 which guarantees that the masses add up to Svr^//?, the instanton action 
per unit (imaginary) time. At higher topological charge k, the parameter space of 
widely separated constituent monopoles is described by kn monopole constituents, k 
of each of the n types of abelian charges (with overall charge neutrality). 

We established in an earlier paper ^Ul that well-separated constituents act as point 
sources for the so-called far field (that is far removed from any of the cores). When 
constituents of opposite charge {n constituents of different type) come together, the 
action density no longer deviates significantly from that of a standard instanton. Its 
scale parameter p is related to the constituent separation d through p'^/P ~ d. Yet, 
the gauge field is vastly different, as is seen from the fact that within the confines 
of the peak there are n locations where two of the eigenvalues of the Polyakov loop 
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coincide [TT| ll2j. thus in some sense varying over the maximal range available (e.g. for 
SU(2) from I2 to — 12), whereas for trivial holonomy only one such location is found. 



Figure 1: In the middle is shown the action density in the plane of the constituents 
at t = for an SU(2) charge 2 caloron with trPoo = 0, in the regime where con- 
stituents are not well-separated. On a scale enhanced by a factor lOvr^ are shown the 
densitities for the two zero-modes, using either periodic (left) or anti-periodic (right) 
boundary conditions in the time direction. This solution is for the so-called "crossed" 
configuration with k = 0.997 and D = 8.753, see Sect. HI for more details. 

On the other hand, when constituents of equal charge come together typically an 
extended core structure is found. This was deduced, in particular for the case of 
charge 2 calorons, from our ability to analytically determine the zero-mode density 
(summed over the two zero-modes implied by the index theorem) in the far field limit, 
neglecting exponential contributions^. In this limit it forms a singular distribution on 
a disc bounded by an ellipse, but approaches two delta functions for well-separated 
like-charge constituents. This zero-mode density only sees constituents of a given 
charge, depending on the boundary condition for the fermions in the imaginary time 
direction, which can be chosen to be a U(l) phase (containing the physically relevant 
choice of anti-periodic boundary conditions for thermal field theory). We will show for 

^This is in some sense equivalent to the high temperature Umit, with constituent masses given by 
Sn'^i'm/P, such that the range of the exponential contributions shrinks inversely proportional with the 
temperature. 



3 



SU(2) that their difference for periodic and anti-periodic boundary conditions coincides 
exactly with the (abehan) charge distribution extracted from the gauge field in the 
same hmit, making contact with an old result due to Hurtubise jTHj for the asymptotic 
behaviour of the monopole Higgs field. 

We found two particular parametrizations within the SU(2) charge 2 moduli space 
that exhibited these extended charge distributions. The first includes as a limit arbi- 
trary charge 2 monopoles. The second of these parametrizations contains as a limiting 
case the axially symmetric configurations constructed for arbitrary charge in Ref. ^3]. 
Deforming away from the axial configuration the two discs overlap. Describing the 
intricate behaviour for the non-abelian core of these configurations in this region of the 
parameter space requires one to find exact solutions, which are presented here. We rely 
on early work of Nahm ^Hj and Panagopoulos JEj for charge 2 magnetic monopoles, 
which is simplified to some extent by our formalism that uses the relation between the 
Fourier transformation of the ADHM construction (as relevant for the finite temper- 
ature case) and the Nahm transformation, a crucial ingredient for our success to find 
explicit solutions . Fig. ^ gives a particular example for the action and zero-mode 
densities of a charge 2 caloron solution. The two dimensional zero-mode basis is chosen 
such that each zero-mode only localizes on one of the constituents of a given charge, 
showing both the zero-modes with periodic and anti-periodic boundary conditions in 
the imaginary time direction. 

This paper is organized as follows. In Sect. El we will outline the construction, 
introduce the Green's function that is computed through the analogy of an impurity 
scattering problem, and summarize the various limits that can be formulated before 
explicitly solving for the Green's function. In Sect. El we present the method that 
allows one to find the exact solution for the Green's function, first for the general case 
and then applied in more detail to that of topological charge 2 calorons. Readers only 
interested in the results could skip Sects. 12.21 and El In Sect. HI we discuss the two 
classes of configurations in the moduli space of the charge 2 calorons and provide plots 
of the various quantities to illustrate the properties of the exact results. In Sect. El we 
discuss the relation between the algebraic tail of the gauge field and the zero-mode 
density. We end with some discussions. An Appendix presents a new result for the 
limiting behaviour of the action density. 

2 Outline of the construction 

There are two steps in the construction of charge k caloron solutions. The first step 
involves finding a U{k) gauge field A^{z), which satisfies the self-duality equation, i.e. 
the Nahm equation on a circle parametrized by z, with z introduced through 
replacing the original SU(n) gauge field Afj,{x) by A^{x) — 27riz(5o^l„. Although not 
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affecting the field strength, this changes the holonomy to exp{—2TTizf3)Voo, revealing 
that z has period (3~^. The index theorem guarantees the existence of k zero-modes 
\&(a;; z) which satisfy the Dirac equation, or in the two-component Weyl form 

L>^*(x; z) = a^D^^ix; z) = a^{d^ + A^{x) - 27nz6o^,ln)'^{x■, z) = 0. (2) 

with (Jfj, = ajj^ = (I2, —ir) {ti are the usual Pauli matrices). We may remove z from 

the gauge field A^j_{x) by transforming the zero-mode to = exp(— 27rzt2;)\l/(a;; z), 

which is at the expense of making the zero-mode only periodic up to a phase factor, 
^zit + = exp(— 27riz/3)\l/2(t, x). In a similar way we could introduce z through 
^{x; z, z) = exp{2niz- x)\l/(a;; z), which replaces in Eq. ^ ^m(^) ^m(^) ~ 27riz^l„, 
where Zq = z. Assuming the k zero-modes z, z) to be orthonormal one has 

Af{z)= [ ¥''\x;z,z)^—'^^''\x;z,z)d'^x, (3) 

oz^ 

or equivalently (demonstrating as well that A does not depend on z) 

Af{z) = J zy—¥^\x; z)<fx, Af{z) = 2ni J ^^''\x; zyxk¥^\x; z)<fx. 

(4) 

We have shown how A^{z) can alternatively be related to a Fourier transformation 
of the ADHM construction of instantons jTTj, that periodically repeat (up to a gauge 
rotation with Poo) in the imaginary time direction, so as to turn an infinite charge 
instanton in to one of finite charge and finite temperature. The derivation of this 
relation will not be repeated here, see Ref. O ^5 for the details. 

The connection to the ADHM construction has been useful to simplify the second 
step in the construction of the caloron solutions, namely how to reconstruct the original 
gauge field when given a solution to the Nahm equation ^H] (which is equivalent to 
the quadratic ADHM constraint), 

^A,{z) + [A^{z),A,{z)\ + \ejM[Mz),Mz)] = 2mY.6{z - fi^m)pL (5) 
az ^ 

For convenience of notation we will henceforth use the classical scale invariance to set 
/3 = 1. The singularities in the Nahm equation appear precisely for those values where 
g-27riz-p^ has one of its eigenvalues equal to 1, at 2; = /ij. This is where some fermion 
field components become massless, i.e. the zero-mode becomes delocalized, whereas 
for generic values of z it is exponentially localized, which has turned out to be a useful 
tool to pinpoint the constituent monopoles. 

One could apply the Nahm transformation again, introducing Afj_{z) — 2iTix^lk 
(with xq = t) and find the n chiral fermion zero-modes in the background of this 
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gauge field. The construction is somewhat comphcated due to the presence of the 
singularities, whose structure is determined from the matrices pm which appear in the 
Nahm equation. Not all pm are independent, e.g. integrating and tracing the Nahm 
equation yields the conditions, Z]m=i Tr pm = 0. Further constraints are implied [Hj 
by the fact that one may introduce (as is most easily seen in relation to the ADHM 
construction) k two-component spinors in the n representation of SU(n), such that 

2TxClP^Cb = (y,St-r-p^, (6) 

where Pm are projections defined through Poo = YJL=i ^'^^^'^'^ Pm {Sm will appear in 
Eq. dHD). 



2.1 The Green's function 

However, with reference to the ADHM construction, there is great benefit in first 
finding the solution for the Green's function, fx{z,z') = g^{z)fx{z,z')g{z'), where 



dz'^ 



+ V{z- X) z!) = 47i\kSiz-z'), (7) 



with 

V{z;x) = 4:Tc'^R'^{z;x) + 27r ^ - /i„)5'„„ = g{Pni)SrngKl^m), 

m 

R,{z- x) = X, - {2m)-'g{z)A,{z)g\z), (8) 

and Sm playing the role of "impurities". This is formulated in the gauge where first 
we transform Ao{z) to a constant (diagonal) matrix, 27ii^o, as is always possible in one 
dimension, and then use 

g{z) = exp(27rz(eo - Xolk)z), Tr^o = (9) 

(when Tr ^0 7^ it is absorbed in a shift of xq) to transform Aq — 27iixolk to zero. This 
is at the expense of introducing periodicity up to a gauge transformation; although 
fx{z, z') is periodic in z and z' with period 1 (for f3 = 1), fx{z, z') no longer is^. 

Given a solution for the Green's function, there are straightforward expressions for 
the gauge field fl] (only involving the Green's function evaluated at the "impurity" 
locations) and the fermion zero- modes fOl ^] • For the zero-mode density this gives 

¥:^{xm\x) = -{2n)-'dlf:\z,z). (10) 



^It is in this respect interesting to note that ^(l) plays the role of the holonomy associated to the 
dual Nahm gauge field A^{z). 
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In this paper we will only have need for the Green's function at z' = z, which formally 
can be expressed as 

f^{z,z) = -An'{{l,,-J^^)-')^^, ^,^^t(i)pexp£^'(^^^0^^^ ^^^w, (11) 

where the (1,2) component on the right-hand side of the first identity is with respect 
to the 2x2 block matrix structure. In particular this leads to a compact expression 
for the action density [3 

Six) ^ -\iiFl^{x) = -\dldl logdet (^e--"«(l2. - T^) , (12) 

which can be shown to be independent of the choice of z. 

The formal expression for Tz can be made explicit by a decomposition into the 
"impurity" contributions Tm at 2; = /i^ and the "propagation" Hm = W^m(/^m+i, yUm) 
between firn and /im+i- For z G (/im,/im+i) this gives 

with 



^--\2.S^ Ij' ''-'^'^'-\m{z)fj-{z))\fj^^{z')fj-{z') 

(14) 

The columns of the two kx k matrices fm{z), defined for z G (/im,yUm+i), form the 2k 
solutions of the homogeneous Green's function equation, 

of which those in f^{z) are exponentially rising and those in f^{z) are exponentially 
falling. This implies [Hj f^{z) exp (±27r|f|(z — /im)lfc) for \x\ oo, in which 
Cm = /m(/^m) cau be arbitrary non-singular (to ensure a complete set of solutions) 
matrices. In Ref. [TU] we put = 1^, but here we find it convenient to leave this 
choice open. With the "impurity" scattering problem solved, constructing the exact 
solutions of the homogeneous Green's function equation is the last step in finding 
analytic expressions for the higher charge calorons. 



2.2 Limiting cases 

Nevertheless, approximate solutions can be derived, either assuming x is far removed 
from any core such that the gauge field has become abelian (which we called the far 
field limit, denoted by a subscript "fF'), or assuming x and the constituents of type m 
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(belonging to the m-th interval) are well separated from all others, but not necessarily 
from each other (which was called the zero-mode limit, denoted by a subscript "zm"). 
This is because we found jTHI that for z G (/im.,/im+i) the k zero-modes \f^"^(x) only 
"see" the constituents of type m. For \im < z < fim+i we have 

{fmif^n.)f;ni^)-' ' ^+/+ (/in^) (;^) i?-'(^), (16) 

up to exponential corrections in the distance to the constituents of type m' ^ m, with 
= Ifc ~ 2S^^i?m,_i(/im), Z+ = Ifc - 2S„^i?m(/im,), (17) 

Rm{z) = ki^mi^) + = R^il^m) + Ki-lil^m) + Sm- 



and 



In the zero-mode limit Z+ and Z^^^ approach 1^ up to algebraic corrections. Neglect- 
ing these contributions as well, e.g. by sending the constituents of type m' ^ m to 
infinity, leads to the so-called monopole limit (denoted by a subscript "mon"), further 
simplifying the expression for the Green's function to 

/^°''(2:, z) = -7rf/(2:,/im+i)f/^^(/im,/im+i)?7m(/im,z)i?~^(z), (19) 

with^ 

u^{z, z') ^ mz)mzr' - (20) 

In turn, from Eq. (fT^ one may derive the far field limit, giving up to exponential 
corrections in the distance to all constituents 

/f(/i„,/x^) = 27rS^i, (21) 



as is relevant for the expression of the gauge field in this limit ^3] . For the zero-mode 
density (Eq. (UHl)) we may use for /i^ < z < /im+i {z strictly different from ^m, l^m+i) 

f^{z,z) = 7iR-^\z). (22) 

The discontinuity in this limit, nRmifim) 7^ SttS^-*^ 7^ irRm-iifim), arises due to the 
zero- mode developing a massless component when z approaches fim- It might seem that 

^Note that Um{z, z') satisfies the Green's function equation with boundary conditions Um{z, z) = 
and £Umiz, z') = -^rUmiz, z') = 27ri?„(z), for z' z. 
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Eq. ()22|). combined with Eq. ()10j). is inconsistent with an exponential decay. However, 
it turns out that 

Vmix) ^ {^Try'Tr {r;^\z)) (23) 

is independent of z in the m-th interval and harmonic everywhere (hence giving vanish- 
ing zero-mode density) except for some singularities in the cores of the constituents of 
type m. It is this feature, and our ability to compute Vm{x) exactly for SU(2) charge 
2 calorons, that allowed us to make statements about the localization of the cores, 
without solving the Green's function exactly. 

In Appendix A we derive the following new result for the monopole limit of the 
action density (Eq. p2|) ). In the limit where x and the constituent locations of type m 
are well separated from all other constituents, for which both the action and zero-mode 
densities are static, we find (with Um defined in Eq. (1^ ) 

5--(f) = -la^aflogdet [f/™(/i„+i,/i„)i?-^(/i^)] , (24) 

up to algebraic corrections. This is a direct generalization for the action density of 
a single BPS monopole, S{x) = —^dfOj log sinh(27rz/m|x|)/|5| (located at the origin 
and with mass Syr^i/m)- We emphasize that this result can be used irrespective of the 
distance between the constituents of the same type m. Eq. (j24j) therefore gives in terms 
of fmi^) ^ closed form expression for the multi- monopole energy density. The same 
holds, when using Eq. (jl9p . for the zero- mode density (Eq. (jiup ). For the caloron it 
involves taking the limit where all constituents of type m' ^ m are sent to infinity, 
which is why this is called the monopole limit. We recall, that in the far field one 
should use [13] (see also App. A) 

n 

S^{x) = -\d]d] Y: logdet (/+(/.^+i)/^(/^m)-^i?-^(/x™)S^^) . (25) 



3 Exact results 



Finding the exact homogeneous solutions of the Green's function equation, Eq. (fTH|l . 
closely follows Nahm's method ^] to construct the dual chiral zero-modes. The main 
advantage of our approach is that we need not worry about boundary conditions, 
as this is solved by the "impurity" scattering formalism ^1]. In the following we 
restrict ourselves to a given interval z G (/imj/Um+i) and work in the gauge where 
Ao{z) - 2Tcixolk = 0. 

Using that A^{z) is self-dual, a consequence of the Nahm equation, one easily shows 
that DIDx = + 4:T^'^R'^{z; x), such that it is natural to consider the equation 

dJ{z) = a^D^Jiz) =(-^- 2nf ■ R{z- x)] i,{z) = 0. (26) 
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It follows that ipi^) would be a homogeneous solution of the Green's function equation, 
albeit that ip{z) is a spinor (with a chirality opposite to that for the zero- modes involved 
in the Nahm transformation, cmp. Eq. (j2|)). We follow Nahm and use the ansatz 
ipi^z) = (I2 + u{x) ■ t)\s) ® v{z), where v{z) is a /c dimensional complex vector, u{x) is 
a unit vector that does not depend on z and \s) is an arbitrary normalized constant 
spinor (as long as it is not annihilated by I2 + u{x) ■ f). It then follows that v{z) = 
{s\ (I2 - u{x) ■ f)'ii){z) satisfies Eq. (fT3|) . 

The unit vector u{x) is found from a complex vector y{x) which squares to 0, 
y{x) ■ y{x) = 0, implying its real and imaginary parts are perpendicular and of equal 
length (7^ as long as y{x) 7^ 0), such that (for ease of notation the x dependence of 
y and u will henceforth be left implicit) 

u = iyxy*/{y-y*) (27) 

is well defined and u x y = —iy, i.e. Re(^), Im(y) and u form an orthogonal set of 
vectors. Using the ansatz for ip{z) and introducing 

Y{z) = -y- R{z;x), U (z) = ~27ru ■ R{z; x) , (28) 

leads to the equations 

Y{z)v{z) = 0, -^Hz) = U{z)v{z) (29) 

for which the first one can only have a solution provided detY{z) = 0. 

It is the great beauty of Nahm's formalism that det Y{z) is a conserved quan- 
tity. That is, -^Aj{z) = —^ejki[Ak{z),A£{z)] (the Nahm equation Eq. (0) restricted 
to an interval and in the gauge where Aq = 0) implies ^ det Y{z) = for any 
choice of x and y on a null-cone in (or rather in CP^ since we may rescale y 
with a non-zero complex factor without changing the equations). The conserved 
quantities are generated by the symmetric traceless monomials, Mi^i^...i^, built from 

[Ai^{z)Ai^{z) ■ ■ ■ Aii,{z)^ with ^ arbitrary, as one readily verifies. For example, 

Ti Ai[z) is constant and defines the center of mass. A natural way to project on 
the traceless symmetric monomials is precisely through introducing y G CP^ on a 
null-cone, forming yi^yi^ ■ ■ -yi^Mi^i^..^^ = yi^yi^ ■ ■■yi.Ti (^Ai^{z)Ai^{z) ■ ■ ■ Ai^{z)y It is 

interesting to note that Xi^^Xi^ ■ ■ ■ Xi^Mi-^i^...ii,\x\~^ is always a spherical harmonic of order 
£, used in Ref. [12] to show through the multipole expansion of Tr (_R~^(z)) that it is 
conserved and harmonic, except for singularities in the core of the constituents. 

Once it is established that detF(z) is independent of z for any choice of a:, we can 
look for its zeros. Using the null-cone parametrization y = (5(1 — C^)) ~t(l + C^); C); 
and the fact that the matrix Y[z) is k dimensional, we obtain a polynomial equation in 
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( of order 2k and hence there are for generic x exactly 2k solutions. It is useful to note 
that these solutions come in complex conjugate pairs, where the symmetry y ^ y* 
implies u —u and C, ^ —1/C,* (giving y* up to a multiple, equivalent to y* in CP^). 

Given a particular zero y, we may conveniently write a vector in the kernel of 
Y{z) as [in] Va{z) = (f){z){sidiY{z))ac for a fixed choice of c, where SidjY^z) is the 
matrix formed by the minors of Y{z). The Nahm equation is easily seen to imply 
j^Y{z) = [U{z),Y{z)] for any choice of y on the null-cone. From this one derives that^ 

^adjf(z) = [U{z), ad] Y{z)]. Substituting Va{z) = (f){z) {ad] Y{z))ac into Eq. (f29|l 
gives 

^(adjF(^)),e = 0W(adjF(^)f/(;.)) . (30) 
az ^ ^ ac 

Using the fact that U{z) = —iujg{z)Aj{z)g''{z) — 2nujXj, we get 

exp( il{z) -27T zu-x) da(z) i{adiY{z),Ujg{z)Aj{z)g^{z)\ 

(f)( = ^ — ^ - ±_ttc_ (3X) 

v/(adjF(.)).e ' dz 2{adiY{z)U 

where the equation for fi{z) is the same for any value of a (it may depend on the value 
of c). This is useful for studying the asymptotic behaviour of the solution. For large 
detY{z) = implies that y-x^O, such that (cmp. Eq. (j2ZI)) u ±x/\x\. We 
may use the symmetry u — > —u to guarantee that there are k zeros y^''^ with the sign 
of M*^^)- X negative, leading to solutions that rise as exp(27rz|a^|). It then follows that 
the k zeros = give ■u(^+'^) = —u^^\ leading to the solutions that decay as 

exp(— 27rz|x|) for large Hence we may put /Jb(-z) = v^\z) and f~f,{z) = v^J^~^''\z). 
Defining 

ni+{z) = -{adiy^'^-R{z;x))^^, m-,{z) = - [ad] f R{z- x)) (32) 
which are algebraic in Aj(z) and xj, we find 

Jli^) = ^.izH^'Hz), Uz) = m-,{z)^^'^^\z), (33) 

where (j){z) contains the exponential dependencies, and thus seemingly all the infor- 
mation about the cores of the constituents. To make this more precise we compute 
R^{z), see Eq. (HHI), using that Eq. ^ implies (with u^''+^^ = -m^^^) fjabi^) = 
=F27rX)^=i ^^'^'^ ■ Rad{z] x)f^{z), finding that the factors (p^^'^z) drop out 

RaM = - E ^^'^ ■ Rae{z-,x)mf,{z){m^{z)-'),,. (34) 

b,e=l 



''Assume first that y is such that det F 7^ 0, in which case adj Y ~ Y ^ det Y and therefore 
^adjf = -Y-^[U,Y]Y-^detY+Y^^-^dctY = [C/, adj y]+f-iTr(y-i^y) = [C/,adjf]. Observe 
that adj Y is analytic in y, such that the result is valid also when y leads to det Y = 0. 
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This proves that R^i^) is purely algebraic in Aj{z) and and Rm{z), 

which determine the far field limit for the zero-mode density and the gauge field. 

One might wonder how, given that the ^jJ^''\z) are of the "wrong" chirality in the 
context of the Nahm transformation, one could use these results to reconstruct the 
gauge field for magnetic monopoles where the relation to the ADHM construction is 
not readily available. For this one observes that the columns of il'^''\z) can be used to 
form a2kx2k matrix w{z). Using that Dx'w{z) = 0, one finds Dl (^w''{z)~^^ = 0. Thus 
the columns of w'^{z)~^ give 2k independent solutions for each interval, from which n 
normalizable solutions '^^p\z; x) should remain after imposing the appropriate bound- 
ary (cq. matching) conditions. These are then used in Nahm's original construction to 
compute the gauge field (cmp. Eq. Q) 

AP;{x) = I ¥P\z;xo,x)^-^¥'^\z;xo,x)dz, (35) 

where '^^p\z;xo,x) = {z)'^^p\z; x). 

There seems to be considerable advantage in using the Green's function (Fourier 
transformed ADHM) method, since it can solve the matching conditions without relying 
on the availability of exact solutions for the normalizable dual zero-modes. To go 
beyond the approximations discussed in Sect. 12.21 and resolve the constituent cores we 
need to solve for fi{z). This can not always be done in closed form, but it is given by 
an explicit integral which can be performed numerically when required. For charge 2 
monopoles Panagopoulos [ini was, however, able to find the exact integral. We can use 
the same ingredients for the caloron case and explicitly solve for the Green's function 
in the case of charge 2 calorons. 



3.1 Analytic expressions for charge 2 

For charge 2 the number of invariants associated to the conserved quantities of the 
Nahm equation is 8, of which TTAj{z) = Aniaj are related to the center of mass for 
the constituents of given magnetic charge, coming from the interval under considera- 
tion. Assuming now that Aj{z) is traceless, 5 invariants remain, given in terms of the 
symmetric traceless tensor 

M,, = -i {TTiMz)A,iz)) - l5,,TTiAliz))) . (36) 

Three of its parameters are associated to the rotation TZ which diagonalizes the 3x3 
matrix, M = 7^diag(ci, C2, Cs)?^*, where 71 is fixed by requiring C2 < ci < C3. The q 
add to zero and can be expressed in terms of the so-called scale (D) and shape (k) 
parameters, 

ci = , C2 = , C3 = D^— — . 37 

1 12 12 12 ^ ^ 
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The Nahm equation for the case of charge 2 can be solved completely in terms of 
Jacobi elliptic functions 1^ , which was summarized in Ref . jTUI 

g{z)Aj{z; a, 7^, h, D, k)g\z) = 2TTiajl2 + ^iDUjbfb iD{z - Zq)) h)nh, (38) 

where /i is a global gauge parameter and 

f ( \ - ^ f ( \ k'snk(^) dnk(z) , r 

/i ^ = TT, f2{z) = h{z) = --, k = VI - k2. 39 

CUkU CUkU CUkU 



with^ snk(2:) = sm{(f{z)), cnk(z) = cos{(f{z)) and dnk(-2) = y 1 — k2sn^(2;) the stan- 
dard Jacobi elliptic functions. This does not yet address the matching of Aj{z) on the 
different intervals, where some difference between the monopole and caloron applica- 
tion appears. For the caloron, apart from the axially symmetric solutions constructed 
in Ref. we found two sets of non-trivial solutions that interpolate between over- 
lapping and well-separated constituents. It is for these classes of solutions that we will 
resolve the cores, when constituents overlap and the non-linearity plays an important 
role. 

The next step in the construction is finding the zeros y of det Y{z). One has 

det Y{z) = y^yj{x,Xj - if '(5^, - {2n)-^Mij), (40) 

where we used y"^ = 0. Substituting a parametrization for this null-cone in CP^, e.g. 
y = (i(l — C^),— 1(1 -|- C^))C)) gives a 4th order polynomial. However, for finding 
the 4 solutions we find it in this case more convenient to first diagonalize the matrix 
XiXj — ^x'^6ij — (27r)~^Mjj = OikXkOjk- Introducing y' = O'^y, the equation for the 
zeros reduces to (y^)^Ai + (2/2)^^2 + (z/3)^A3 = 0, which in addition to the null-cone 
condition, (y'l)^ + (^2)^ + iVs)^ = 0) is now easily solved by 



^i-)'= (^^A2-A3,(-l)'^+V^3-Ai,VA2-Aij, y(-+^)'=(y('^)y, a = 1, 2, (41) 

where we fixed (for generic x) the diagonalizing rotation O by ordering Ai < A3 < A2. 
Using Eq. (jTfjl we find for u' = O^u 



u 



V V A2 — Ai VA2 — Ai / 



One easily checks that ■u'^^'^^ and -u*^^'^) give rise to respectively the exponentially rising 
and falling solutions. 

It is also instructive to give the explicit expressions for the matrices m^{z) in 
Eq. (jn21) (choosing c = 2). We note that for a 2 x 2 matrix adj F = (TrF)l2 - F, 



(p{z) implicitly defined by the elliptic integral of the first kind z = 1/v 1 — sin^ 6 dO 
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and without loss in generality^ we take zq = 0, a = 0, 71 = l^,, D = 1 and /i = I2 in 
Eq. dng), such that 

^> An [Arrx . yW - yi'^ Mz) At^x ■ - yf' h{z) ) ' 



We checked that Eq. ^,V{x) = (27r)-iTr(i?+(2)+i?-(;z))-\ evaluated using Eq. 
is independent of z and agrees with the result derived in Ref. [TU] . 

We next address solving Eq. (|3ip. which for charge 2 can be written as^ 

dfljz) ^ MjjUjyjd"'' - 27ri{x ■ y)uiA^''{z) 
dz 277 (f ■ - iyiAfiz) 

For the same choice of parameters in Eq. (pH|l as above, zo = 0,a = 0,7^ = l3,-D = l 
and h = I2, this gives (with a = 1 and c = 2) 

^/^(^) _g^^^ -.s ^l/l(^)-m2/2(^) . . 

Although the dependence on x is complicated, the integral over z turns out to be 
manageable (as was observed before in the context of charge 2 monopoles, although we 
here choose not to express the solution in terms of theta functions JEI)- To solve the 
equation we first rewrite the right-hand side of Eq. ()45|1 using the fact that uxy = —iy 
(cmp. Eq. (j2ID), 

'"1/1(2;) - iu2f2{z) {usyi + iy2)fiiz) - i{u^y2 - Wi)f2iz) 



yifiiz) - iy2f2iz) 2/3 (2/1/1(2;) - 22/2/2(2;)) 

/i(^)/2(^)2/| + 422/i2/2(kr , ^3 
2/3 (167r2(f ■ y)2 - ylfi{z)) 2/3 ' 

In the last identity we used that y"^ = 0, /f (2;) — /K-z) = 1 — k^ = (k')^ and the fact that 
detF(2;) = implies (2/1/1(2;) - iy2f2{z)){yifi{z) +^2/2/2(2;)) = IQti'^{x ■ yf -ylfl{z). 
With -^hiz) = fi{z)f2{z) we can now integrate Eq. (jl^ . 

AW = 2..(2 ■y)^ + i log ( ft ■ i^^ ) + ,!!S!lM(^Mi/(.). 

2/3 \^T^x ■ y - fsizjys J 2it{x ■ y) Ay^ 

I{z) ^ Yi^{f^\z),n)-Yi^{l,n) + \zl n ^ (46) 

1/3 

^We may change zq, a, ft,, 7?. and by resp. translations, (gauge) rotations, and suitable rescalings. 
'^Using that ^{Ai,Aj}yiUj = ^12^^ [Ai A j)yiUj = -Mijy^Ujl2, since m • y = 0. 
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up to an irrelevant constant, where nk(s, n) is the elhptic integral of the third kind^, 

nk(s,n) ^ / , (47) 

•^0 (1 - nt^)J{l-W){l-t^) 



We now combine these ingredients to give in terms of y the exact form for the 
homogeneous solution of the Green's function equation, 

Va{z) = 0(^)(adj Y{z)U = exp (/i(^) - 27rzu ■ x) (43) 

v/(adjr(2))i2 

or putting in all the relevant expressions 

Va{z) = expfz-[27r(fxf)3 + /i^^^^(k|+nk(/3-i(^),n)-nfc(l,n))l') X 
\ 2/3 L 87r(a; ■ y)\z\ ^ 

(A \-^( f(\ ( . / x\3 /47rf-y- (-l)'^?/3/3(z)y 

Substituting y = yC"^ = O^y^^'^', with y^^'^' as defined in Eq. dHJ gives f^^{z) = v^J'^z) 
and f~i,{z) = v^^^^{z), and thereby the Green's function, once we specify the parame- 
ters involved in the solutions to the Nahm equation. 

4 Action and zero-mode density plots 

The discontinuities in Aj{z) at z = fii and z = ^2 implied by the Nahm equation, 
Eq. (0), impose constraints which are (like the quadratic ADHM constraint) in general 
difficult to solve. Work is in progress to describe the full parameter space for SU(2) 
and charge 2, but in Ref. JU] we did find two non-trivial parametrizations for which 
we illustrate here in a number of figures how the full caloron solutions look like, using 
the exact Green's function as constructed in the previous section. Taking advantage of 
the possibility to work with arbitrary arithmetic precision the programme Maple has 
been used for these calculations. The configurations are formulated in terms of the 
two intervals z G [/xi,/^2] and z G [/i2, 1 + /^i], each associated with two constituent 
monopoles of equal magnetic charge, but opposite in sign from one interval to the 
next. Apart from a shift and (gauge) orientation, the configuration is described by a 

^More commonly the elliptic integral of the third kind is defined as n(n; Lp, k) = Ilk (sin n). Note 
that I{z) can alternatively be written as j/'^'^'' {l—t^/n)^^{t^ — \<L^)^'^{t^ — l)^^dt, from which it follows 
that {-J{z) = ((1 - nfl{z))\h{z)h{z)\)-' ±h{z) = sign(z)(l - n/|(z))-i, using /|(z) - 1 = /|(z), 

_d_ 

dz " 



fi{z) - k2 ^ fUz) and ^J,{z) = h{z)h{z) 
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shape (k) and scale (D) parameter (see Sect. I3.1|) . for simplicity assumed to be the 
same on both intervals. We also take all constituents to be of equal mass, 1/1^2 = 5 
(/i2 = — yWi = j), most relevant for the confined phase with trPoo = 0. 

The two periodic and two anti-periodic chiral fermion zero-modes each have support 
on oppositely charged constituents (see Fig. and in the far field limit it was found 
that the zero-mode density (summed over the two zero-modes implied by the index 
theorem) is described by a disc singularity, bounded by an ellipse with semi-major 
axis D/An and eccentricity (1 — k')/(l + k') (where k' = -y/l — k^). This revealed 
that the core of a cluster of like-charged constituents is in general extended, unless the 
individual constituents are well separated. The far field only describes the (algebraic) 
abelian component of the gauge field and to resolve the structure of the core we need 
to determine the full non-abelian structure. 



y 




Figure 2: An example to illustrate the location of the disc singularities (light and 
dark shaded according to magnetic charge) for a rectangular (left) and crossed (right) 
configuration. The latter is shown at a = — 7r/2, for d = 71/32 and 6 = tt/A (k = 0.962, 
D = 3.894 and ip = — 7r/4). The curves indicate the would-be constituent locations 
at fixed d and 6, varying a from — tt (to which the arrows point) to 0. For a = —it 
and the discs collapse to lines (k = 1) with no singularities remaining, except at the 
endpoints. 

For the first parametrization (called "rectangular") the two discs are parallel and 
separated in height by a distance d. The configuration is charaterized by the two 
extremal points along the major axis of the disc, also called would-be constituent loca- 
tions^, 

y^^ = {0,k{-ird , i-iyiAn)-'D) , 2nd = DU\D), (50) 

^In their immediate neighbourhood the action density is maximal; they are the constituent locations 
in the point-like limit, k — s- 1. 
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up to an overall shift and orientation (the definition of f2{z), which involves k, can 
be found in Eq. ()39|) ). A typical example is shown in Fig. |21 (left). Apart from k and 
D, the parameters that enter Eq. ()38|1 for the m-th interval are a = (0, i(— 0), 
71 = I3, h = I2 and zq = ^{1 + (—1)^) . In all cases discussed here we have also put 
1^0 = 0. One verifies that the discontinuities of Aj{z) at z = fj,m are given by 2nip{^ 
with appropriately chosen (a (cmp. Eq. (jH))), as discussed in detail in Ref. [TUj . 



Figure 3: Shown is the energy density (middle) for fi2 = 0.25, k = 0.570 and D = 6.915, 
in the monopole limit d ^ 00 and in the x-z plane through one of the discs for the 
rectangular configuration, see Fig. |21 (left). On a scale enhanced by a factor Att"^ are 
shown the densities for the two monopole zero- modes (left and right). 

For the second parametrization (called "crossed") the two discs are coplanar and 
intersect, see Fig. |2] (right) for a typical example. Their relative orientations can vary 
between perpendicular and coinciding (for which k is forced to 1). Here we choose for 
the parameters in Eq. ()38j) h and TZ to be non-trivial (isospin) rotations around the 
y-axis with angles (— l)'"^, resp. (— 1)™V9, and a = (0,0, — §(— l)™'(icosa), whereas zq 
and ^0 are as in the rectangular case. The would-be constituent locations are now given 
(up to an overall shift and orientation) by 

yO) = ((-l)i(47r)-iDsin^, 0, (-l)"+^(47r)-iD cos^ - i(-l)'"dcosa) , (51) 

where ±(p conveniently gives the orientation of each of the two discs with respect to 
the 2-axis. The angle a originates from the definition of (a, which through Eq. © 
determines the discontinuity of Aj{z). To ensure the proper matching, the following 
three equations need to be satisfied jTU] 

L)sin(^±y?)[/3(iD)±/i(iD)] = 87rd(l±sina), Df2{\D) + 87idsma = 0. (52) 
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which determine y?, k and D for given a, 9 and d. Fig. |21 illustrates that for these 
crossed configurations the two discs always overlap, unlike for the rectangular case. 
This formed an important motivation for the present study, so as to determine in how 
far the overlapping discs would affect the behaviour in the core. 



Figure 4: In the middle is shown the action density in the plane of the constituents 
at t = for an SU(2) charge 2 caloron with trPoo = in the crossed configuration of 
Fig. 121 hence with k = 0.962 and D = 3.894. On a scale enhanced by a factor IGvr^ are 
shown the densities for the two zero-modes, using either periodic (left) or anti-periodic 
(right) boundary conditions in the time direction. 

We will illustrate, using the exact formalism, how the action and zero-mode densities 
for these solutions behave. For the rectangular case we are mostly interested in the 
monopole limit, d oo. For finite d, the density separates in two contributions, 
where each in the limit d cxd is exactly the density for a charge 2 monopole with 
the same values of k and D. Note that in the limit d ^ oo the matching conditions 
for Aj{z) "decouple" the intervals, turning into pole conditions at z = fii^2, as is 
appropriate for multi-monopoles [201 ■ We do recover from this the known results for 
the charge 2 monopoles like the doughnut structure for k = 0, corresponding 
to two superimposed monopoles. The interest in the monopole limit comes from the 
fact that the zero-mode densities for multi-monopoles had not been studied in detail 
before In Fig- El we give the densities for k = 0.570 and D = 6.915, which is 

intermediate between the doughnut and well-separated monopole configurations. On 
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the other hand, for d small (compared to /? = 1) the configuration will look like two 
non-dissociated calorons, and in particular is no longer static. When D remains much 
bigger than d, forcing k — > 1, these behave as two well-separated charge 1 instantons. 
Otherwise, when D is comparable to d, one finds overlapping instantons 

An example for the crossed configuration with k = 0.997 and D = 8.753 was already 
shown in Fig. ^ In this case D is large enough for the like-charge constituents to be 
separated, as is particularly clear from the zero-modes, which are essentially no longer 
overlapping. But two nearest neighbour (oppositely) charged constituents still show 
appreciable overlap. The distance between these nearest neighbours is \D\/2/t{ = 
0.985. As this is comparable to /3 = 1 we would expect the configuration to depend 
on time. Indeed, at the maxima of the action density its value of 1.18 x IGvr^ at t = 
is reduced by almost 50% at t = 0.5. At the center of mass, where the action density 
is much lower, there is still a time dependence. However, far from all cores the field 
becomes static^''. Increasing D further (which will push k closer to 1) the configuration 
quickly turns into well separated spherically symmetric static BPS monopoles. 

More interesting is to consider the case with smaller D, like D = 3.894 and k = 
0.962, for which the disc singularities were illustrated in Fig. |21 The corresponding 
densities are shown in Fig. 0] We see that the constituents are now so close that they 
form a doughnut, but we stress this is different from the static monopole doughnut, 
which has k = 0. Since the oppositely charged constituents now are as close as 0.438, 
which is considerably smaller than the time extent, the solution will have a strong time 
dependence. When D is decreased even further, it will turn into a charge 2 instanton 
localized in both space and time. 



5 Higgs field asymptotics 

In this section we make some comments on the far field limit of the gauge field. As we 
discussed before, the gauge field far removed from any core becomes abelian (as well 
as static). The abelian subgroup is the one that leaves the holonomy invariant, in the 
periodic gauge equivalent to leaving the constant asymptotic value of the adjoint Higgs 
field Aq invariant. For definiteness, let us consider the case of SU(2) with k = 2, f3 = 1 
and Voo = exp(27ria; ■ f) (i.e. /i2 = —fii = \uj\). Up to exponential corrections we have 

A^{x) = 2muj ■ f - ^iuj ■ T^{x), (53) 

where we can express $(a;) in terms of the far field limit of the Green's function at the 
impurities ^1] 

$(x) = n-' \l - n-'Tr (/^^(/xa, 1-^2)82)] 9,Tr (/^^(/ia, /X2)pQ • (54) 



^°For any SU(n) and topological charge k the far field limit is static, provided all constituent 
monopoles have a non- vanishing mass 



19 



Using the twistor description of magnetic monopoles Hurtubise [12] was able to 
explicitly compute the asymptotic Higgs field for the SU(2) magnetic monopole long 
ago. The function he found for this algebraic tail amazingly agrees exactly with Vm{x), 
Eq. (j2Sl), which was introduced to describe the caloron zero- mode density (m denotes 
the interval and hence the type of constituents to which the corresponding zero-modes 
would localize JHl)- As mentioned before, from our multi-caloron results one can 
recover the multi-monopole results by sending the constituent monopoles with the 
"unwanted" magnetic charge to infinity, cmp. Eq. 

Although this tends to be cumbersome to show, in the far field can be written 
as $i(x) — $2(3;), where ^m{x) is the contribution coming from the type m constituent 
monopoles and the difference in sign is due to the sign change in the magnetic charge. 
This is simply because the field is abelian in the far field and linear superposition 
preserves the self-duality. Hence, ^m{x) = 27rVm(^), such that for the SU(2) caloron 

^(f) = 27rVi(f) - 27rV2(f). (55) 

We checked that this result indeed holds for the solutions discussed in Sect. El even when 
the two types of monopole structures are not well-separated. This relation trivially 
holds for the axially symmetric solutions that were introduced in Ref. jT3], where $(x) 
was explicitly shown to factorize in a sum of point charge contributions, compatible 
with what was found in Ref. ^lOj for the zero-mode densities. Therefore, in the far field 
limit (i.e. for the algebraic part) the singularity structure in the zero-mode density 
agrees exactly with the abelian charge distribution, as given by df^{x). Such a relation 
is at the heart of using chiral fermion zero-modes as a filter to isolate the underlying 
topological lumps from rough lattice Monte Carlo configurations 0. 

6 Discussions 

In this paper we have analyzed the higher charge caloron solutions and showed how 
to obtain exact results by suitably combining techniques developed in the context 
of the Nahm transformation and the ADHM formalism. The aim of these studies 
has been to establish that SU(?7.) caloron solutions of charge k can be described in 
terms of kn monopole constituents, and that these can be viewed as independent 
constituents. A natural way to get an ensemble would be to consider approximate 
superpositions of k charge 1 calorons, but this would lead to an unwanted memory 
effect, with constituents remembering from which caloron they originated |14,. Our 
studies, within the context of self-dual configurations, have shown nevertheless that 
the constituents have an independent identity, with the only requirement that the 
net magnetic and electric charge of the configuration vanishes (each of the n types of 
constituents should occur with the same number). A recent lattice study 1^, using 
the technique of over- improvement [^H]; fully confirms this picture. 
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It is therefore reasonable to consider the constituents as the independent building 
blocks for constructing an ensemble of monopole constituents, something that was 
not questioned in Ref. 0, but like for the instanton liquid |22j forms an essential 
assumption in a semi-classical study. Clearly the expectation is that semi-classical 
methods no longer work in the confined regime, at least for the part of the parameter 
space that corresponds to well-separated constituents, that is typically associated to 
instantons with a large scale parameter. It is not unlikely that the density of these 
constituents at low temperatures is so high that they form a coherent background and 
as such will no longer easily be recognized as lumps. With high quark densities leading 
to deconfinement, it may perhaps be that a high constituent monopole density will 
lead to confinement |2HI; although for now we have to leave this as a speculation. 

Instantons that overlap get deformed and depending on the relative gauge orien- 
tation tend to "repel", i.e. inspecting the action density distribution they do not get 
closer than a certain distance j^lj. When deconstructing instantons in monopole con- 
stituents, interestingly only like-charge constituents will show this effect, manifesting 
itself through the extended core structure. For unlike charges, from the point of view 
of the abelian field, the configuration behaves as with linear superposition. If as a 
consequence of this all abelian charge is annihilated, it disappears through forming a 
small instanton (localized in space and time), which in the limit of zero size describes 
the boundary of the moduli-space. The interaction between constituents of opposite 
duality is more complicated [SI 123 • 

In conclusion, calorons with non-trivial holonomy have revealed a rich structure, 
incorporating traditional instanton physics, but allowing for gauge fields that inherit 
some essential features associated to a confining background not present in the tradi- 
tional formulations. The fact that the underlying constituents are monopoles opens the 
way to describe the confining aspects of the theory in terms of these degrees of freedom. 
Much work remains to be done when it comes to understanding the dynamics, but we 
hope to have convinced the reader that a consistent picture is developing that holds 
considerable promise for the future. 
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Appendix A 



In this appendix we derive the zero-mode hmit for the action density, which assumes 
that the distance of x and the constituents of type m to all constituents of type m' ^ m 
is large, but where x and the constituent locations of type m may otherwise be arbitrary. 

As in Ref. [Iflj we take z = Hm + for computing JF^ and use that we can write 
det(l2fe-J^;.„) = det(l2fc-Lir), where K = Fm-iQm-i ■ ■ ■eig^{l)Fnen- ■ ■Qm+2Fm+l 
and L = 6m+i-Fm0m, with 



27ri?+(/im) -27ri?„(/im) / " V27ri?+_i(/im) -27ri?.^„i(/i„ 

fmif^m+l) fmif^rn) ^ \ ^^^^ 

/m(/^m,+l)/m(/^»n)"V 

We note the K has no remaining dependence on the constituent locations of type m. 
Writing LK = LK + LK, with 



we find det(l2fe - LK) = det{K) det{k-^ -L- LKR-^). 
We next use 

-ilKvi -K~lK^\ ~^y_, ( 



(57) 



and note that in the zero-mode limit K^^, K^^Kj^^^, K_^K^^ and (K^^) are 

exponentially small (cmp. Ref. jTU], App. A), such that 

det (ze-™«(l2fc - LK)) = det(e-2""°;r++) det(L++). (59) 
With the definition of L we now find 

L++ = jRm+lifJ'm+l) (-Rm+l(/^m+l) + 'S'm+l) Um (^R^_i{flm) + Sm) , (60) 

where 

f^m = 2m+l/m(A''m+l)/m(A''ni) if^rn) — Z^_^ifj^(flrn+l) fmif^m) R-m if^rn)^^, 

2:^ = ^k±{R~{f^m) + Sm)~^R^_^{fin,), = Ifc ± i?^(^„) (i?+ _i (/i„) + S„) (61) 

and Um contains all contributions due to the constituent locations of type m, up to 
exponential corrections in the distance of these, and of x, to the other constituents. 
Hence log det {ie~'^'^^° {l2k — LK)) splits into the sum of two contributions, logdet(f/m) 
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and logdet[i(i?+_i(/irri) + Sm)e^'''^"K++R^^\^{jj,m+i){Rm+i{l^m+i) + Sm+i)], where the 
last term only depends on the constituent locations of type m! ^ m whose contribution 
will decay inversely proportional to the fourth power of their distance. Allowing for 
algebraic decay (or in the monopole limit, sending all constituents of type m' ^ m io 
infinity) such that in addition 2^^+! = = '^k, one thus finds Eq. (j^ . 

A simple way to derive the result for the far field limit in Eq. is by noting 
that in this case all fmif^m+i) fmif^rn)~^ are exponentially small and Fm can be ap- 
proximated by diag(F++, 0), with F^~^ = fmifJ^m+i) fmif^rn)^^ ■ This therefore acts as a 
projection on the ++ component and is thus seen to lead to det{ie~'^^^°{l2k — LK)) = 
det(e-2-«o^t(i)i7;++e++...i7^++e++). Using the fact that HOj 9++ = ^R-\nm)^m 
gives the required result. 
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